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I I Abstract 

Q We calculate the Fermi surface of electrons in hole-doped, extended t-J models on a square 

^ lattice in a regime where no long-range antiferromagnetic order is present, and no symmetries 

are broken. Using the "spinon-dopon" formalism of Ribeiro and Wen, we show that short-range 
antiferromagnetic correlations lead to a reconstruction of the Fermi surface into hole pockets which 
I are not necessarily centered at the antiferromagnetic Brillouin zone boundary. The Brillouin zone 

area enclosed by the Fermi surface is proportional to the density of dopants away from half-filling, in 

O 

> 

phase transition from the FL* state to the Fermi liquid state with long-range antiferromagnetic 
O order. 

O 
(N 



C/3 



contrast to the conventional Luttinger theorem which counts the total electron density. This state 
realizes a "fractionalized Fermi liquid" (FL*), which has been proposed as a possible ground-state of 
the underdoped cuprates; we note connections to recent experiments. We also discuss the quantum 
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I. INTRODUCTION 



The nature of electronic Fermi surfaces in strongly correlated metals, in particular under- 
doped cuprates, has been the subject of intensive debate for many years. Recent observations 
of pocket-like Fermi surfaces in quantum oscillation experiments^"^ as well as new angle re- 
solved photo emission measurements^ have triggered a renewed theoretical interest in this 
matter 

One possible, well known route to a Fermi surface reconstruction is the onset of spin- 
density wave (SDW) order, which breaks a large Fermi surface into small electron- and 
hole-pockets centered at the magnetic Brillouin zone boundary^°'^^. In fact, many of the un- 
resolved theoretical problems in strongly correlated electron materials, from heavy-Fermion 
compounds to high-Tg cuprates, are related to the fate of electronic excitations close to anti- 
ferromagnetic quantum critical points^^. It has been been argued, however, that the critical 
point between a metal with a large Fermi surface and an antiferromagnetic metal with small 
Fermi pockets may be replaced by a new intermediate phase, the so called fractionalized 
Fermi liquid (FL*)^^'^^, which exhibits small pockets similar to the antiferromagnetic metal, 
but breaks no symmetries: summaries of these arguments, and of previous theoretical work, 
can be found in two recent reviews. ^^'^^ 

The simplest picture of the FL* phase appears in the context of Kondo lattice models 
coupling a lattice of localized / moments and a conduction band of itinerant c electrons. 
There are two important energy scales to consider: the Kondo exchange Jk between the 
/ moments and the c electrons, and Heisenberg exchange Jh between the / moments. If 
Jk ^ Jh, then / moments are "Kondo-screened" by the conduction electrons, leading to 
a Fermi liquid ground state with Fermi surfaces enclosing the traditional Luttinger volume 
which counts the density of both the / and c electrons; the only memory of the localized 
nature of the underlying / electrons is the that electronic quasiparticles near the Fermi 
surface have an effective mass which is much larger than the bare electron mass, and so this 
phase is often referred to as a 'heavy' Fermi liquid. However, in the opposite parameter 
regime Jh 3> Jk other phases can appear. The most natural possibility is the appearance 
of magnetic order of the / moments, but let us assume the /-/ couplings are sufficiently 
frustrated so that this does not happen. Then the / moments may form a spin liquid, which 
does not break any symmetries of the lattice Hamiltonian. The formation of the spin liquid 
also quenches the Kondo effect, and so the effective value of Jk does not renormalize to 
infinity as it does in the single impurity Kondo model^^'^^'^^. The c electrons are now only 
weakly coupled to the / spin liquid, and so the c electrons form a "small" Fermi surface 
which encloses a volume controlled only by the density of c electrons, which violates the 
traditional Luttinger count. This is the FL* metal. 

This paper describes a FL* state in a sm^'Ze-band model appropriate for the cuprate 
superconductors. Previous studies^^'^^ realized such a state by initially fractionalizing the 
electron into a neutral S — 1/2 spinon, and a spinless "holon" carrying electromagnetic 
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charge. The spinons eventually became the excitations of a 'background' spin liquid, analo- 
gous the spin liquid of the / electrons above. And the holons eventually captured a spinon to 
reconstitute as electron-like particles which occupied the states inside a small Fermi surface. 
Because of this somewhat intricate sequence of transformations, the description of the FL* 
state was only achieved in a semi-phenomcnological manner. 

Here we will provide a more direct and quantitative description of the FL* state in a 
single-band model. The key step will be a rewriting of the single band degrees of freedom 
in a manner which mimics those of the Kondo lattice. Such a formulation is provided by 
the representation of Ribeiro and Wen^° in which the electron fractionalizes into a neutral 
spinon and a "dopon" which has the same quantum numbers as the electron. 

The rest of the paper is organized as follows. In section II we introduce the extended t-J 
model in the representation of Ribeiro and Wen^°, which is ideally suited for our purposes. 
Section III deals with our approach to construct FL* ground states and presents results for 
the shape and position of the electronic Fermi surface. Section IV describes the Fermi surface 
evolution from the FL* state to the Fermi hquid state with long-range antiferromagnetic 
order, along with a discussion of the quantum-critical properties. We summarize our results, 
and note connections to recent experiments in Section V. 



II. MODEL 

In the following we want to study ground-states of extended t-J Hamiltonians on the 
square lattice 

H=-]^Yl^ij (clcja + h.c) + \Y1 -^^^ ■ ~ l^i^j) ' (2-1) 

ij ij ^ ^ 

where c|g. (Cj^) denotes the Gutzwiller projected creation (annihilation) operator of electrons 
with spin a on lattice site i, Si — c^ia<^apc^p is the electron spin operator and rii — cl^c^^^ 
the electron number operator (here and in the following we sum over repeated spin indices). 
We are interested in describing possible ground states slightly below half filling n — 1 — x, 
where the density of doped holes is small a; <^ 1 but large enough to destroy any long-range 
magnetic order. In addition these ground states should not break any lattice symmetries. In 
particular we want to show that strong short-range antiferromagnetic correlations already 
lead to a reconstructed Fermi surface consisting of small hole pockets, the area of which 
is proportional to the dopant density x, instead of 1 — a; as for conventional Fermi liquids. 
Such ground states realize a fractionalized Fermi liquid^^'^*^. 

Our starting point is the spinon-dopon formulation of the t-J model developed by Ribeiro 
and Wen^°. In this representation the elementary excitations are spinons, which carry spin- 
1/2 but no charge, and dopons, which carry spin-1/2 and charge. Accordingly, Ribeiro 
and Wen introduce two degrees of freedom per lattice site, a 'locahzed' spin-1/2 as well as 
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TABLE I: Single site basis-state correspondence: t-J v.s. spinon-dopon. 



t-J 


spinon-dopon 


It). 








|0)i 


(im-i;t).)/%/2 


unphys. 


triplet-states 


unphys. 


doubly occupied dopon 



a fermionic spin-1/2 degree of freedom - the dopon - representing a doped charge carrier. 
A physical hole corresponds to a singlet of a lattice spin and a dopon. The correspon- 
dence between single-site basis states is shown in Tab. I. Following this approach, the t-J 
Hamiltonian in Equ. (2.1) takes the form^° 



^ = ^ E ^^^(S, • S, - 1/4) V{1 - dldJi^ - d%d,^)V 

+1 E y^[i4^,a - lidl^apdj^) ■ (S, + S,) + dld,^S, ■ Sj 

ij 

+ i (Aa^apdj^) ■ (S. X S,) + h.c] V-^iJ2 dUia ■ (2-2) 



Here V — Ylji^ ~ ^jt^it^j-l-^jJ.) denotes the Gutzwiller projector for the fermionic spin-1/2 
operators d| and di that create or annihilate a dopon on lattice site i, and we added a 
chemical potential /i for the dopons. Note again that the spins Sj on each lattice site j are 
independent, localized spin-1/2 degrees of freedom and are not associated with the spin of 
the dopons. This representation of the t-J model is faithful in the sense that the Hamiltonian 
doesn't couple the physical singlet- and the unphysical triplet states in the enlarged Hilbert 
space that is spanned by the spin- and the dopon degree of freedom^". A projection to the 
physical Hilbert space is thus not necessary. 

In terms of the spin- and dopon operators the Gutzwiller projected electron operators 
take the form 

cl = ((1/2 + aS])d,., - SJd,,) V , (2.3) 

where S'^ denotes the spin raising (lowering) operator S'^ {S') for a (X). Prom Equ. (2.3) 
one can easily show that total density of electrons is given by 

E 4c,, = P(l - 5: n"" 1 - Y: d]A„ . (2.4) 

i.e. the density of doped charge carriers equals the density of dopons x — X]o.(^L^io-), as 
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expected. 

As will be explained in more detail below, our main assumption is that the localized lattice 
spins form a Z2 spin liquid with bosonic spinon excitations. This is reasonably justified in 
the doping regime close to the antiferromagnetically ordered phase, where the interaction 
between the lattice spins is frustrated by the motion of dopons. The bosonic nature of 
the spinous prohibits a hybridization of spinous with fermionic dopons and gives rise to an 
electronic Fermi surface, the volume of which is determined by the density dopons x alone, 
as long as no pairing instabilities occur. This is in contrast to the conventional Luttinger 
theorem, which states that in a metal without broken symmetries the 'volume' enclosed by 
the Fermi surface should be proportional to the total density of electrons 1—x. It has been 
shown earlier, however, that topological excitations associated with the emergent gauge field 
of a spin liquid have to be included in the Luttinger count giving rise to a FL* phase with 
small pocket Fermi surfaces the total volume of which is the same as in an antiferromagnetic 
metal. In the present formalism, this modified Luttinger theorem of a Fermi surface of size x 
can be easily proved by applying the usual many-body formalism to the system of interacting 
spinons and dopons described by Eq. (2.2) (and more explicitly in Eq. (3.3) below); we need 
only assume that the final state is adiabatically connected to a state of weakly interacting 
spinons and dopons, and then the standard proof leads here to the novel Luttinger count of 

X. 

In the presence of strong local AF correlations the most important couplings between the 
dopons and the localized spins arc the two interaction terms in the second line of Equ. (2.2). 
The ~ Si ■ Sj term leads to a strong suppression of nearest neighbor hopping of dopons in 
a locally AF ordered background, whereas the ~ (Sj + Sj) term is responsible for scattering 
of dopons with momentum transfer close to q = (11,71). In the following we are going to 
neglect the ~ Sj x Sj term due to the expected strong local coUinear AF correlations. Also, 
we use a mean- field decouphng of the Heisenberg exchange term in the first line of (2.2), 
i.e. J — > (1 — x^J and drop the Gutzwiller projectors, which is safe in the low doping limit 
X < 1. 

III. FL* AND ELECTRON FERMI SURFACE 

As mentioned above, the main prerequisite in order to get a fractionalized Fermi liquid is 
that the localized spins form a spin liquid. Within our model (2.2) spin liquid ground states 
can be conveniently described using a Schwinger-Boson representation for the lattice spins, 
i.e. we write ^ 

Si = ^bl^t^a^b^^ , (3.1) 

which requires the constraint J2a ^la^ia = 1 to hold on every lattice site. Note that there is 
an emergent U{1) gauge structure associated with the redundancy of the Schwinger Boson 
representation under local phase transformations bj^ — > bjaexp{i(l)j). The exchange terms 
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can be expressed in terms of Schwinger Bosons using the identity 

1 \ 5 ■ 



(3.2) 



Inserting these expressions in the Hamiltonian (2.2) and using the approximations mentioned 
above, we obtain the Hamiltonian 
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^ocAab]^ e^sbi^bjs 



(3.3) 



+ 2 E f [2 4ad,a - dldjp {bl^b,^ + b^^b^^) + h.c] + (A blb,^ - J , 



where A is the Lagrange multipher that enforces the Schwinger Boson constraint on average. 
In this representation a spin liquid can be conveniently described by employing the mean- 
field decoupling 

Q^j = li^apblb}^) . (3.4) 

By construction this mean-field decoupling preserves the SU{2) invariance since Qij is a 
singlet expectation value. After a Fourier transformation we obtain the euclidean mean- 
field action (we use the shorthand notation /c = (a;, k)) 



k,a 



-Ep<9p>^p-k if^n + A 



~ XI dk'+g-ka ^k ^k^kq d^>^> + ^ Qt+qJqQk 



q,k,k' 



kq 



•7^ (^k' + ^k'+q+k) [BkiBq^dk+k'+q^-dk'i + Bl^B*^dk'idk'+k+(j^] , (3.5) 



q,k,k' 



where we have introduced the bosonic Nambu spinor 



Bk 



* 



(3.6) 



and 



as well as 



* k'kq — 



^(tk' + ^k'+q-k)5(T,t J2p Qpitp+k'-k + J^p-k'-q) 

Z^p Qp(^p+k'-k + ^p-k'-q) K^k' + ik'+q-k)'5a,i 



Ck — ^k + 2 ^2 Qq^pQk+q-p " 1^ , 

qp 



Saa' , (3.7) 



(3.8) 
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where is the usual tt't" tight-binding dispersion on the square lattice with non-zero hopping 
amplitudes up to third nearest neighbors 



tk = ^ ( cos kx + cos kyj +t' 2 cos cos ky + t" ( cos 2kx + cos 2kyj . (3.9) 

Note that in contrast to Jk and Qk we have absorbed a factor 1/2 in the definition of ik- 
Moreover, ioun and iQ„ denote fermionic and bosonic Matsubara frequencies, respectively. 
Because the bosonic spinon modes are gapped in the spin-liquid phase we can safely integrate 
them out and obtain an effective action for the dopon fields d. Expanding to second order 
in the bosonic propagator we get 

S^^Ud, d] = Sold, d] + Trlog/3^;o"' - Tr6?o$ 

- ^Tr Qo^Qo^ - Tr g^DG^D - IV GoDQ^D , (3.10) 

where Tr denotes the trace with respect to momentum-, Matsubara- and Nambu indices. 
Furthermore we have defined 



(^0 ^)kq - 



(3.11) 



and 



^kq — C^fc'+g-fco-Vk'kq^fcV (3-12) 

k' 

Dkq — J y^(^k' + ^k'+q+k)^fc'.j,'^fc'+fc+gt x (3.13) 
Dkq 



4 
k' 



4 ^(^k' + tki+ci+l!.)dk'+k+q^dk'i. (Tx , (3-14) 
k' 

where denotes the respective Pauli matrix in Nambu space. The effective dopon action 
in Equ. (3.10) describes the hopping of doped charge carriers in a locally AF-ordered back- 
ground as well as the residual interactions between dopons due to the exchange of a spinon 
pair. Note, however, that the Schwinger boson mean-field theory presented here cannot 
be used to describe a conventional Fermi liquid state with a large Fermi surface at large 
doping, where the Luttinger volume is determined by the total density of electrons 1 — x. 
This phase can be described using Schwinger Fermions instead of Bosons^°. In this case a 
hybridization between dopons and Fermionic spinous can lead to a ground-state with a large 
Fermi surface. 
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A. Gaussian dopon action 



The linear contribution ~ Qo to the effective dopon action (3.10) basically descends from 
the term in the Hamiltonian (2.2) which couples the dopons to Sj ■ S-,- and thus strongly 
renormalizes the bare dopon dispersion. Performing the Matsubara summation and the 
trace over the Nambu indices of the Tr Qq^ term we get the Gaussian action 



k,a k' 

+'IVlog/36;o'' + -^q-kQ^Qk . 

kq 



d 



ka 



(3.15) 



This expression can be simplified using the self-consistency conditions for the Lagrange 
multiplier A as well as for the mean-field Qk- We make one further approximation here, 
however, and determine A and Qk at the Gaussian level not fully self-consistent, but only 
with respect to S^p = Trlog/3^Q"^, i.e. we neglect the back-action of the dopons on the 
spinons. The approximate self-consistency equations thus read 



1 = 



dX 



E 



A 



x^^p-k 

k 



Qk- 



X]q "^q-kQq 



2Eu 



(3.16) 
(3.17) 



where is the free energy associated with S^p. Here and in the following denotes the 
spinon dispersion relation, which is given by 



(3.18) 



Inserting these expressions back into (3.15) we get 

^MF ^ /^^^dkai-it^n + ^k] dka + Tr log /J^g"^ + COUSt. 
k,a 

with the Gaussian dopon dispersion 

Ck — — 2 ^ QqQpip+k-q — ■ 

qp 



(3.19) 



(3.20) 



Note the different sign of the second term compared to the bare dopon dispersion in 
Equ. (3.8). 
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B. Mean-field ansatz for a Z2-FL* 

In the following wc use the simplest mean field ansatz for the Qi/s. We take a zero- 
flux state with Qij = Q on nearest neighbor bonds. The invariant gauge group^^ (IGG), 
i.e. the group of gauge transformations that leaves this ansatz invariant, is U{1) due to 
the bipartite nature of the square lattice. Indeed, we can choose the gauge transformation 
bja- bjcr exp{i(f)) on sublattice A and bja — >■ bj^ exp{—i(f)) on sublattice B without changing 
the ansatz. However, since gapless U (1) gauge fluctuations are hard to control, we break the 
U{1) gauge group down to Z2 by including frustration in the form of a small next-nearest 
neighbor exchange interaction J' as well as a corresponding singlet bond-amplitude Q'. The 
excitations of the emergent Z2 gauge fleld are gapped visons which shouldn't play a big role 
in our subsequent analysis as long as their gap is sufficiently large, thus we neglect them in 
the following. The Fourier transformation of our ansatz Qij then takes the form 

Qk = i2 [(5(sin kx + sin ky) + Q'{sm{kx + ky) -\- s,m.{ky — kx))] ■ (3-21) 

Note that Qji = —Qij and thus Q-k = —Qk- The corresponding spinon dispersion relation 
(3.18) is given by 

E], = ^^X^ - ^JQ{suikx + sm.ky) + 2J'Q' coskx suiky]"^ , (3.22) 

where J and J' are the nearest- and next-nearest- neighbor exchange couplings whereas Q and 
Q' are the corresponding singlet amplitudes on nearest- and next-nearest neighbor bonds. 
For this choice of Qk the convolutions in all the expressions for the effective action can be 
evaluated straightforwardly. In particular, the Gaussian dopon dispersion from Eq. (3.20) 
takes the form 

= - A\Q\'^) [cos kx + cos ky) 

t" 

-\-t\l — 4|(5'|^) coskxcosky -\- — (cos2A;^ -|- cos2ky) — ji . (3.23) 

Note that the dispersion is invariant under Z2 gauge transformations Qij — )■ —Qij. The 
singlet amplitudes Q and Q' can take values between Q, Q' G [O, 1/a/2] , where Q = Xj \/2 if 
nearest neighbor spins form a singlet. One can clearly see that the nearest neighbor hopping 
amplitude vanishes for perfect classical local AF correlations (Q=l/2) and it changes sign 
for Q > 1/2. It is important to emphasize, however, that the Gaussian dopon dispersion 
(3.23) is strongly renormalized by the residual interaction. 
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FIG. 1: Dopon-dopon interactions induced by the exchange of a spinon pair. 

C. Residual interactions, effective dopon action at quartic order 

Here we analyze the interactions between dopons that are induced by the exchange of a 
spinon pair. The quadratic terms Qq in the effective dopon action (3.10) are given by 

SinI = - I] Tr2 [l{go)k^k<i{Go)q^qk + {go)kDkgigo)gD<ik + {go)kDk,{go)<iD,k\ . (3.24) 

kq 

The first term gives rise to non-spinflip interactions, which we denote by V^^ic'' whereas the 
second and third terms describe interactions where the dopon spins are flipped (denoted by 
^kk')- Both interactions are shown schematically in Fig. 1. 

In the following we perform a self-consistent Hartree-Fock analysis of this induced re- 
tarded interaction. Self-consistency is necessary because the interactions are strong and the 
shape as well as the position of the Fermi surface is strongly affected by interaction induced 
fluctuations^^. The Hartree-type interactions are accounted for already to a large extent in 
the Gaussian dopon dispersion (3.23). In fact, the Hartree diagrams would correspond to a 
self-energy correction of the bosonic spinon propagator and are not expected to change the 
results qualitatively. For this reason we restrict ourselves to the two Fock-type diagrams 
shown in Fig. 2. Since we expect the Fermi-liquid character of the dopons to prevail, we 
use a dominant pole approximation and neglect the incoherent part of the dopon Green's 
function. In this approximation, the dopon Green's function in the diagrams in Fig. 2 takes 
the form 

G^iKiu)^—^^ , (3.25) 

where the quasiparicle residue as well as the dopon dispersion ,^k are calculated self- 
consistently. We will justify this approximation a posteriori by checking that the quasipar- 
ticle weight Zk is reasonably large. The dopon self-energy corresponding to the diagrams in 
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FIG. 2: Fock-type self-energy diagrams for the dopons. 



Fig. 2 thus takes the form 

E(k,ia;„) = 



k+q 



T/(l) , T/(2) 

^qkk ^ ^qkk 



(3.26) 



+ ^^9) Okk'q + ^'kk'q 



where akk'q and 6kk'q are momentum dependent factors given by 



(3.27) 



Okk'q = 3(tk + tk+q)^ — 2|((5*t)k'-k + (Q*^)k'+k+q| 
&kk'q = (tk+tk+q)'[3A2+ ((g*J)kKg*J)k'+q + C.C.)] 

+ A(tk + tk+q) [{Q**J)k' + (Q** J)k'+q] [{Q*t)k'-k + (g*t)k'+k+q] + C.C. 
+ {Q**J)k'{Q**J)k'+q[{Q*t)k'-k+ (Q*i)k'+k+q] +C.C. 

+2A'|(g*i)k'-k + (Q*i)k'+k+qr . (3.28) 

Here the asterisk denotes convolutions, i.e. {Q*t)i^ = XlpQp^k-p- After performing the 
Matsubara sums and analytic continuation ia;„ (jj + i6 we get (at T = 0) 



Sj?(k,a;) 



5]^k+q j 
q,k' I 



-©(-a+q) 



-E'k'-E'k'+qakk'q — &kk'q 
+q a; — £'k' — -Bk'+q — ^k+q + 

(_Ek' — Ck+q + l^)akk'q + ^kk'q/ -E'k' 



El^^ - {E^' - a+q + + ^'^)' 

(-£^k'+q + Ck+q ~ '^)akk'q + ^kk'q/ -E'k'+q 
El, - (-Ek'+q + Ck+q -UJ- 



(3.29) 



The imaginary part of the retarded dopon self energy finally takes the form 

&kk'q Okk'i 



ImSi?(k,a;) 



TT 

16 



k'q 



k+q 



2£Jk'-E'k'+q 



p) [e(a+q)5((^ - i?k' - i?k'+q - a+q) 



+e(-ek+q)5(a; + Ek' + £;k'+q - Ck+q)] 



(3.30) 
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Note that IioErIuj) = for — 2A < a; < 2A, where A denotes the spinon gap. 

In all subsequent calculations we do not determine the Lagrange multiplier A (which fixes 
the Schwinger-Boson constraint) self-consistently, but use it to fix the value of the spinon gap 
A. Moreover, we use the nearest- and next-nearest neighbor singlet amplitudes Q and Q' as 
free parameters. Our calculational procedure works as follows: first, we calculate lmSi^(k, u) 
numerically using the adaptive Monte-Carlo integration algorithm MlSER^^, which is based 
on a recursive stratified sampling method. Since the computational effort increases consid- 
erably with increasing accuracy, we set the bound of the relative error estimate to be smaller 
than 6%, which is arguably a relatively large value, but sufficient for our purpose. In order 
to perform the Monte-Carlo integration we smoothen the singularities of the delta-functions 
as well as the step-functions by replacing the delta-functions by Lorentzians with a FWHM 
of 0.01 and the step functions by Fermi distributions at an effective inverse temperature 
P = 200. 

The second step is to evaluate the real part of the self energy by a Kramers-Kronig 
transform and determine the dopon dispersion by finding the maximum of the dopon 
spectral function 



An \ -ImEij(k,a;) 

A{k,uj) = . (3.31) 

^ [-a; + + ReS«(k,a;)]2 + [ImSR(k,a;)]2 



Here, denotes the Gaussian dopon dispersion from Equ. (3.23). The quasiparticle residue 
is obtained via 



(3.32) 



Finally, the self consistency loop is performed by plugging ^k and Zk back into Equ. (3.30) 
and repeating the steps above until convergence is achieved. 

D. Relation between the electron- and the dopon Fermi-surface 

Using Equ. (2.3) the electron momentum distribution can be expressed in terms of the 
lattice spin- and dopon operators as 

cLck. = const. + e^'^-(^-^^) [-(1/4 + S, • S.O^d,, + {dy^^d,^) ■ (S, + S,)/2 

(3.33) 

where we implicitly sum over repeated spin indices and again neglect the Sj x Sj term. 
Using the Schwinger-Boson representation (3.1) for the lattice spins the electron momentum 
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FIG. 3: (Color online) Self-consistent dopon dispersion (left) and dopon quasiparticle residue Zk 
(right) in the Z2-FL* phase as a function of kx and ky in the upper right quadrant of the Brillouin 
zone. The thick black contour marks the position of the dopon Fermi surface, which coincides 
with the electron Fermi surface. The dashed line indicates the magnetic Brillouin zone boundary. 
Parameter values for this plot are: Q = 0.4, A = 0.025, /i = 0.083 and the rest as in Tab. II. 

distribution is given by 

(cLck.) = const. + 1 J]e^k.(R,-R,)^ _ ^t^^^^ + l{e^,blb]^){e,sb,,b^s)d]A. 

= — -(c^^ko-'^-ko-) + smooth function of k . (3.34) 

The second hne follows because the last two terms give rise to convolutions of the do- 
pon momentum distribution with spinon correlators, where the Fermi surface singularity is 
smoothened out. We thus conclude that the electron Fermi surface coincides with the dopon 
Fermi surface. Moreover, the value of the electron quasiparticle residue at the Fermi surface 
is one half times the dopon quasiparticle residue Z^. 

From Equ. (3.34) it is also clear that the spinon-dopon approach to the t-J model cannot 
give rise to electron pockets at the antinodal regions, as observed in recent experiments. 
Instead the Fermi surfaces are always hole-like. 

E. Results 

The following results were obtained using standard values for the bare hopping ampli- 
tudes, shown in Tab. II. The nearest- neighbor hopping amplitude defines our energy scale 
and has been set to unity. The next nearest neighbor exchange interaction J' as well as 
the corresponding singlet amplitude Q' were chosen to be relatively small compared to the 
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FIG. 4: (Color online) As in Fig. 3 but with parameters: Q = 0.54, A = 0.025, = 0.18. 



TABLE II: parameter values 



0.25 


Q 


variable 




variable 


0.05 


Q' 


0.1 


A 


variable 



t" 0.1 



nearest neighbor values, as their only purpose is to break the IGG form U{1) down to Z2. 

Our results for the self- consistent dopon dispersion and quasiparticle residue as 
a function of and ky in the upper right quadrant of the Brillouin zone are shown in 
Figs. 3, 4 and 5. These results are at a finite dopon density > 0, although we note 
that as well as look qualitatively similar in the case where the dopon density is 
going to zero. The position of the dispersion minimum depends on the strength of the 
local antiferromagnetic correlations, which is parametrized by the singlet bond-amplitude 
Q. For Q = 1/2 the hole pockets are aligned with the magnetic Brillouin zone boundary and 
centered at q (7r/2,7r/2). For weaker correlations Q < 0.5 the hole-pockets are shifted to 
the outer side of the magnetic Brillouin zone boundary towards q = (tt, tt), whereas stronger 
local AF-correlations Q > 0.5 give rise to hole-pockets centered on the inner side of the 
magnetic Brillouin zone boundary (see Fig. 6). 

The spinon gap A does not influence the position of the pockets, but it changes their 
shape slightly. The smaller A is, the more elliptical are the hole-pockets. This is illustrated 
in Fig. 7. We note, however, that the ellipticity of the hole-pockets depends more strongly 
on the precise value of the bare hopping parameters t' and t" as on the size of the spinon 
gap. In fact, smaller t' and t" give rise to more elliptical pockets, similar to the standard 
SDW theory for antiferromagnetic metals. 

The effective mass of the dopons turns out to be enhanced compared to the bare electron 
band mass as well. For the two dispersions shown in Fig. 7, the arithmetic mean of the 
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FIG. 5: (Color online) As in Fig. 3 but with parameters: Q = 0.5, A = 0.01, = 0.182. 



effective masses at the dispersion minimum along the two principal axes is rfieg ^ 2.5 in 
natural units (i.e. m = 1 corresponds to the band mass of the nearest neighbor tight- 
binding dispersion). Again the effective mass depends on the precise value of the bare 
hopping parameters t' and t". For highly elliptical pockets the effective mass can reach 
values on the order of nies ~ 10 along the fiat direction. 

Within our approximation scheme the quasiparticle residue does not drop sharply on 
the outer half of the Fermi surface, as expected from phenomenological models^^. Only for 
relatively high dopon fillings, as in Fig. 3, an asymmetry of between the inner- and outer 
side of the hole-poclcet appears. 

IV. PHASE TRANSITION BETWEEN A Z2-FL* AND AN AF-METAL 

The Schwinger-Boson description (3.1) is well suited to study the quantum phase tran- 
sition between the Z2-FL* described above and a metal with long-range antiferromagnetic 
order. Indeed, AF-ordering corresponds to a condensation of Schwinger Bosons at the points 
where the spinon-gap closes. For Q' = the ordering wave vector is commensurate and 
the spinon dispersion (3.22) has two degenerate minima at the momenta q = ±K with 
K = (7r/2,7r/2). A condensate of the two Schwinger-Boson flavors at these respective mo- 
menta, i.e. {bq^) = A/m^5q,K and {bqi) = i/^^q-K, corresponds to an AF-ordered state 
with staggered magnetization mg in x-direction. 

Within our path integral formulation in Equ. (3.5) such a Schwinger-Boson condensate 
can be straightforwardly introduced by shifting the Nambu fields 



and keeping only the quadratic terms in the shifted field. Note that the mean field ansatz 




(4.1) 
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FIG. 6: (Color online) Evolution of the Fermi surface in the Z2-FL* phase as a function of the 
nearest neighbor singlet amplitude Q G [0,l/\/2]. Shown is the upper right quadrant of the 
Brillouin zone. Black solid line: Q = 0.54, A = 0.025, /i = 0.18, red dotted line: Q = 0.5, A = 
0.01, /X = 0.182, blue dashed line: Q = 0.4, A = 0.025,^ = 0.083, other parameters as in Tab. II. 
The hole-pockets move to the inner side of the magnetic Brillouin zone boundary (indicated by the 
dashed line) as the strength of local antiferromagnetic correlations, parametrized by Q, increases. 

Qk in Equ. (3.21) also acquires a contribution from the condensate 

Qk = ^(VK-4,-K) + QL°^ (4.2) 

where Q^f'' = i2Q{suikx + smky) describes the strong nearest-neighbor correlations on top 
of the uniform long-range correlations induced by the condensate. The mean-field action in 
the AF-ordered phase has the same form as Equ. (3.5) with three differences: first, is 
given by Equ. (4.2). Second, the bare dopon dispersion in Equ. (3.8) is replaced by 

= tk(l - m,) - mlit^ - tk-.) + 2 Qt^* QIcVp - ^' (4-3) 

qp 

and most importantly, the condensate gives rise to an additional term to the action (3.5) 
which describes the scattering of dopons with momentum transfer q = tt and which takes 
the form 

^af//5 = — ^ X^(^k + tk+,r)((i<..„,k+7rtf^t.„,k; + h-c) (4.4) 

Now we can perform the same analysis as in Sec. Ill by integrating out the bosonic modes 
and performing a Hartree-Fock analysis of the effective quartic dopon action. At this level 
of approximation the off-diagonal elements of the self-energy in spin-space vanish identically 
and the diagonal elements have exactly the same structure as in Sec. III. The effective dopon 
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FIG. 7: (Color online) Evolution of the dopon dispersion in ttie Z2-FL* phase as a function of 
the spinon gap A for Q = 0.54 and n = 0.15. Left: A = 0.025, right: A = 0.005, other parameters 
as in Tab. II. The Fermi energy is below the dopon band in both cases. With decreasing spinon 
gap A the dispersion around the minima becomes more elliptical. 



action in the AF-ordered phase including the self-energy corrections is thus given by 



-iUn + ^k+T 



(tk + ^k+7r) 



(^k + ik+7 
-iUn + ^k 



d, 



'a;„,k+7rt 



d, 



cj,i,k4. 



(4.5) 



where ^k and Zk again denote the self-consistently determined dopon dispersion and quasi- 
particle residue, calculated in the same manner as in Sec. Ill using Eqs. (4.2) and (4.3). 
Diagonalizing (4.5) we obtain two dopon bands with dispersions 



<.,(k) = ^±i!i±l7te 



k+TT - ^k)^ + m2ZkZk+7r (tk + ^k+Tr)^ • 



(4.6) 



The dopon Fermi surface is determined by w±(k) = 0. Slightly beyond the AF critical 
point, where the condensate density is small (m^ ^ 1), the self energy contributions to 
the dopon Green's function are basically the same as in the Z2-FL* phase with a vanishing 
spinon gap and the dispersion ^k has a minimum close to q = (7r/2,7r/2), as in Sec. III. 
The dopon Fermi surface in the AF-ordered phase thus again takes the form of pockets 
close to q = (7r/2,7r/2). However, the pockets are symmetric with respect to the magnetic 
Brillouin zone (BZ) boundary due to the presence of long-range antiferromagnetic order. 
The shape of the hole-pocket again depends on the strength of the short-range correlations, 
parametrized by Q in Equ. (4.2). Possible Fermi-surfaces for = 0.05 are shown in Fig. 8. 
For these plots we used the same ^k and Zk as in the FL* phase with a vanishingly small 
spinon gap and Q' = 0, which is justified for ^ 1, as argued above. Note that for 
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FIG. 8: (Color online) Possible Fermi surface shapes in the antiferromagnetically ordered phase 
for nis <C 1. Shown is the upper right quadrant of the Brillouin zone. The dashed line marks the 
magnetic Brillouin zone boundary. See text for a discussion. 

(5 = 1/2 the dispersion is almost symmetric with respect to the magnetic Brillouin zone 
boundary and thus the two dopon bands in the AF-ordered phase are almost degenerate for 
rris <^ 1. In this case we get two concentric Fermi pockets, shown as solid blue and dashed 
red lines in the left plot of Fig. 8. For larger m^, the red Fermi pocket shrinks to zero, 
and we eventually obtain the familiar single hole pocket centered on q = (7r/2,7r/2) of the 
AF ordered phase. The right plot shows a Fermi surface for Q = 0.54, in which case the 
single pocket on the inner side of the magnetic BZ in the FL* phase is "symmetrized" at 
the magnetic BZ boundary. 

Note that in the AF-ordered phase the electron Fermi surface is related to the dopon 
Fermi surface not by the same Equ. (3.34) as in the FL* phase. The spinon condensate 
gives rise to additional contributions ~ rris, which only change the electron quasiparticle 
residue, however. The shape of the electron Fermi surface still coincides with the dopon 
Fermi surface. 

The nature of the quantum critical point between the AF-ordered and FL* phases can be 
addressed by methods similar to earlier work^^"^^. The magnetic fluctuations are described 
by the spinor Schwinger bosons, and their critical theory is the 0(4) Wilson-Fisher fixed 
point. We now have to check if this critical point is destabilized by the dopon Fermi surfaces. 
Because the dopons don't carry emergent gauge charges, they couple rather weakly to the 
critical spin fluctuations^^; the influence of this coupling can be analyzed perturbatively, 
and as in previous work^^'^^ it is found to be irrelevant. So the critical theory remains that 
of the deconfined 0(4) variety.^'' 
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V. CONCLUSIONS 



This paper has presented a microscopic construction of a FL* phase in a single-band t-J 
model on the square lattice, and described its evolution towards the onset of antiferromag- 
netic order. This was achieved by writing the t-J model in a Kondo-like formulation using 
the spinon-dopon formalism. The FL* phase had a "background" spin liquid, which was 
the Z2 spin liquid with bosonic spinon excitations. This spinon was then coupled to mobile 
carriers (the "dopons" ) which had the same quantum numbers of the electron. Our effective 
Hamiltonian for the spinons and dopons was an exact, in principle, representation of the 
t-J model. However, we only analyzed this effective Hamiltonian in a relatively straightfor- 
ward self-consistent one- loop approximation. But there is an obstacle to extending such an 
analysis to higher orders and accuracy by using more powerful computational methods. 

Despite its uncontrolled nature, our analysis yielded physically sensible results for the 
electron spectrum, which resemble aspects of the experimental observations. The key feature 
was the presence of a small hole pocket centered near but not at the magnetic Brillouin 
zone boundary. This pocket enclosed a volume determined by ,x, the density of doped 
carriers alone. The quasiparticle residue was anisotropic around the Fermi surface, but our 
approximation did not yield the strong variation found in earlier phenomenological models.^® 

We note a recent independent study^*^ to describing the under-doped cuprates as a 
"Luttinger-volume violating Fermi liquid" (LvvFL) with a spin liquid of fermionic spinons. 
The LvvFL state is qualitatively the same as the FL* state. 

On the experimental front, there are a number of recent indications that FL*-like a model 
of pocket Fermi surfaces without antiferromagnetic long-range order may be appropriate 
for the pseudogap region of the hole-doped cuprates. The angle dependence of quantum 
oscillations in YBa2Cu306.59 is consistent^ with the absence of spin-density wave ordering. 
NMR measurements^^ on YBa2Cu302; have so far not seen antiferromagnetic order at fields 
as high as 30 Tesla, but do report evidence of charge-ordering. Such a charge ordering can be 
superposed on our FL* analysis in a straightforward manner; as long as the charge ordering 
wavevector does not connect the hole pockets of the FL* state, there will be little change 
in the Fermi surface configuration. And we have already noted photoemission evidence for 
pocket Fermi surfaces.^ 
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